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FLUTTER OF CURVED AND FLAT SANDWICH PANELS 


SUBJECTED TO SUPERSONIC FLOW 
By John A. McElman 


SUMMARY 


A small- deflect ion theory is used in conjunction with a Galerkin technique 
to determine the flutter characteristics of simply supported sandwich panels, 
both flat and slightly curved. The sandwich cover plates are assumed to be of 
the same material and thickness and the core material is assumed isotropic. 

The flow Is assumed supersonic over one surface, and static aerodynamic strip 
theory is used to describe the aerodynamic loading. A two-mode solution is used 
to illustrate the effects of curvature, core transverse shear stiffness, and 
midplane loads. It is found that the effects of curvature can be separated from 
the core stiffness effects. 


INTRODUCTION 


Considerations of dynamic aeroelastic instability, or flutter, of skin 
panels are important in the design of supersonic and hypersonic vehicles. A 
substantial body of literature treats various aspects of the panel flutter prob- 
lem. (See refs. 1 to 5 and references cited therein.) Previous work has focused 
attention primarily on homogeneous isotropic panels. The sandwich panel, on the 
other hand, is being utilized increasingly in the design of aerospace vehicles. 
Under most circumstances, sandwich construction provides higher strength-weight 
ratios in a structure than a homogeneous panel and may provide suitable protec- 
tion against micrometeoroid penetration. The primary distinction between sand- 
wich and homogeneous panels is the greater Importance of transverse shear flexi- 
bility in the sandwich panel (ref. 6). The purpose of the present investigation 
is to provide some estimate of the effect of transverse shear stiffness on the 
flutter boundaries for flat and curved sandwich panels. In this investigation 
both face sheets of the sandwich panel are assumed to be of the same material 
and the same thickness and the core material is assumed isotropic. The flow is 
assumed supersonic over one surface, and linearized static aerodynamic strip 
theory is used to describe the aerodynamic loading. 

An equilibrium equation governing the lateral deflection of a curved sand- 
wich panel is derived in reference 7 on the basis of small-deflection theory. 

This equation with modifications to account for the lateral inertia and aero- 
dynamic loadings is used in this analysis. Effects of midplane loads, trans- 
verse shear stiffness, and panel curvature are evaluated. 



SYMBOLS 


A = R x - 2v 2 

a length of panel, in. 


B 

b 

C 

D Q 

D s 


= + V 2 R . 
y 

width of panel, in. 

2t c ,E„a^ 

curvature parameter, — — - — 

D s r 2 ^ 

transverse shear stiffness of isotropic sandwich panel, lb/in. 
flexural stiffness of isotropic sandwich panel, 

E sts 3 ^ E st s h 2 . _ 

— r + ru xn-lb 

6(l - p 2 ) 2(1 - p 2 ) 


E g Young's modulus for faces of sandwich panel, lb/sq in. 

h depth of sandwich panel measured between middle surfaces of 

faces, in. 

j,m,n,p,f integers 

K jn = J 2 (j + !) 2 + v 2 (2j 2 + 2j + l) + v^ 
k 2 = 7a\o 2 / D s 

Ljn = (2j 2 + 2j + l) + 2v 2 

l lateral aerodynamic loading, lb/sq in. 

M Mach number 

M m p aerodynamic coupling parameter 

N x ,Ny midplane force intensities, lb/in. 
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core stiffness parameter. 


dynamic pressure, , lb/sq in. 


radius of curvature, in. 
time, sec 

thickness of sandwich face plates, in. 


velocity of airflow, in. /sec 
lateral deflection of panel, in. 
amplitude of mn mode, in. 


Cartesian coordinates 


mass per unit area of panel, lb-sec^/ii 


dynamic-pressure parameter. 


critical value of A 


Poisson’s ratio 


mass density of air, lb- sec 2 /in .* 1 
circular frequency, radians/sec 
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A comma followed by a subscript denotes differentiation of the primary 
symbol with respect to the subscript. 


STATEMENT OF THE PROBLEM 


The configuration to be analyzed 
consists of a simply supported curved 
sandwich panel of uniform radius of cur- 
vature r, mounted in a rigid wall with 
air flowing over the top surface at Mach 
number M. (See fig. 1.) The face 
plates are of the same material and are 
of constant thickness t s . The core has 

a uniform thickness h - t s , where h 
is the distance between the middle sur- 
faces of the face plates. The panel has 
a length a and a width b and is sub- 
jected to constant midplane force inten- 
sities N x and Ny (positive in com- 
pression) . The in-plane shear force 
intensity is assumed to be zero. 



The equilibrium equation and appro— Figure 1.- Panel and coordinate system, 

priate boundary conditions are as fol- 
lows (see refs. 1 and 7 ) * 


(la) 

(lb) 

(lc) 

where 

V 2 + 

dx 2 Sy 2 

dx dx 2 c)y 2 dy^ 


.A ♦(!-!£ V* * ♦ * & ♦ * ^ ♦ 7 0 - W,t)' 


w(x,0,t) = w(x/b,t) = w(0,y,t) = w(a,y,t) = 0 


v (x,°,t) = v = w (0,y,t) = w (a,y,t) = 0 
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The inverse operator is defined "by V ^ (v^w) = v^(v \r) = w. The 

remaining terms in equation (la) are defined as follows: 


E c 


n 


H 


B c 


v(x,y,t) 

l(x,y,t) 


Young* s modulus for face plates 

transverse shear stiffness of panel 

flexural stiffness of isotropic sandwich panel 

mass of panel per unit area 
lateral deflection of panel 

lateral loading per unit area due to aerodynamic pressures 


For static strip theory, the lateral loading is given by the simple Ackeret 
value 


= - 2a ^ ( 2 ) 

P dx 


where q is the dynamic pressure pU 2 /2, and (3 = VM 2 - 1. The development of 
equation (la) is based on the assumptions that the deflections are small com- 
pared to the total panel thickness, the total thickness is small compared to 
the radius of curvature, and the sandwich core carries no direct stress. 


ANALYSIS 


A solution which satisfies the boundary conditions for simply supported 
edges is assumed as follows: 


w(x,y,t) 


-L 

-I 


W. 


mn 


m=l 



(3) 


where w^ is the amplitude of the mn mode and cd is the circular frequency. 

The frequency cjd is, in general, complex; however, attention is directed pri- 
marily to real values for which the motion is harmonic. Substituting equa- 
tion (3) into equation (la) and nondimensionalizing by multiplying by 

results in the following equation: 


a' 
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c 




+ 


CQ 

(m 2 + v 2 ) 




m 


(Rx + QR x v 2 - 2v 2 


+ Qv^ 
f 

♦X 


+ Qp)m 2 - (b + QBv 2 + Qv^)j> sin^EE^ 
^(Qm 2 + Qv 2 + l) % cos(^] = 0 


W 


where 


C = 


2tc.Eca^ 


■'S' 1 - I S C 

D g r 2 jA 


Q = 


D s ^ 

D Q a 2 


*x 


N x a 2 

it 2 D 0 


R .5£ 
y *%>« 


A = 


2qa3 

dX 


2 _ 7aXx> 2 


k = 


v = 


D s 

na 


B = ^ + v 2 Ry - v 1 *- 

Ji 


The quantities C and Q are dimensionless curvature and core stiffness 
parameters, respectively; R x and Ry are midplane loading parameters; A is a 
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dynamic-pressure parameter; k 2 is a frequency parameter; and v is a length- 
width- ratio parameter. 

Applying the Galerkin technique by multiplying equation (4-) by sin 

a 

and integrating yields the following set of equations for the coefficients w mn : 



- (R x - 2v 2 + QR x v 2 + QB + Qv^m 2 

+ Qv 2 + l)w pn = 0 

(m = 1 , 2 , . . f) (5) 


where 


^mp ~ 


4mp 

p2_ m 2 


LO 


where p+m Is odd 


where p+m is even 


It ean he noted that equations (4) and (5) are uncoupled in the cross-stream 
direction; n is only a parameter of the solution. 


A general two-mode solution is pursued corresponding to any two consecutive 
modes j and j + 1, and a characteristic equation is obtained by setting the 
determinant of the coefficients of equations ( 5 ) equal to zero. By use of the 
definition of flutter discussed in reference 1, the characteristic equation is 
solved for A, and A is maximized with respect to the frequency parameter B 
to give a critical value A cr , which is a measure of the flutter speed. This 

procedure results in the following closed-form expression for A cr : 


'cr 


_ (23 + i) 2 « 4 
83(3 + l) 



K jn^ + L jn 

K jn« 2 + V* + 1 


C [ij a ( j + l) 2 v 2 + ( 2j 2 + 2J + l) 


K 


■jn 
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where 


K jn = j 2 (j + !) 2 + v 2 ( 2 J 2 + 2j + l) + v^ 
L j n = (2j 2 + 2j + l) + 2v 2 


A solution containing two consecutive inodes, j and j + 1, is assumed since, 
in the region of interest, the minimum flutter speed occurs with the coalescence 
of the frequencies of consecutive modes. More accurate solutions can he obtained 
by including more modes or by solving the uncoupled differential equation exactly 
with n as a parameter. The relative accuracy of these methods as applied to 
a flat homogeneous isotropic plate is shown in reference 1. 


If j = 1, equation (6) gives the following two-term approximation for A cr 
corresponding to the first two modes in the streamwise direction and any one 
mode in the cross-stream direction: 


'cr 


= 2zl 
1 6 




K mQ + L ln 
K lnQ 2 + L ln Q + 1 


c(8v 2 + > ] 


K- 


In 


( 7 ) 


where 


K -) n = b + 5v 2 + 


L ln = 5 + 2v 2 


LIMITING CASES 


If the core thickness of the sandwich panel h - t s is allowed to approach 
zero as Dq approaches infinity and if the thickness and flexural stiffness of 
the resulting homogeneous isotropic panel are denoted by t ^ and D, respec- 
tively, then 

h = t s 

H - 2t S 

Q = 0 
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*sH 

12(1 - n 2 ) 


_ 12aHl ~ H 2 ) 

and equation ( 7 ) reduces to 

( 8 ) 

where A = - 2v 2 . Equation (8) corresponds to equation 1.21 in reference 2 

(where m and r in equation 1.21 are equal to 1 and 2, respectively) for the 
flutter of a curved homogeneous panel. Further, by letting the panel radius r 
approach infinity, equation (8) becomes 

*cr ’ |£<5 - A) (9) 

which corresponds to equation (l6) in reference 1 for the flutter of a flat 
homogeneous panel. 


= 9t±L - A) + ~ H 2 )(5v^ + 8v 2 ) 

16 jAt^r 2 ^ + 5v 2 + v 1 *-) 2 


RESULTS AND DISCUSSION 


Flat Sandwich Panel 


The critical dynamic -pres sure parameter 
can be obtained from equation (6) by letting 


A cr f°r a flat sandwich panel 
C = 0. Equation (6) then becomes 


t 

\ 


x cr 


_ (2j + l) 2 it 4 

83(J + l) 




^jn^ + -^jn 
K jn^ 2 + L jnQ + 1 


( 10 ) 


The value of A cr in equation (l0) is a function of j, v, and Q. For 

a given value of Q and a moderate value of a/b, the flutter boundary (A cr as 
a function of R x j is a minimum for j = n = 1, for the range of con- 

sidered herein. 
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To illustrate the effects of the core stiffness parameter Q, a value of 
a/b of 1 is chosen. For j = v = 1, equation (10) becomes 


A 


cr 


= 9^/p + 10ft + 7 \ 

1 6 y ioq 2 + 7 Q + 1/ 


(ID 


The flutter boundaries ^A cr as a function of R x j for various values of Q 

are presented in figure 2. The boundary for Q = 0 is the boundary for a flat 
homogeneous panel of thickness 2t s . The exact solution for a flat homogeneous 
panel is also shown in figure 2. (See ref. 1.) The flutter boundary is very 
sensitive to the value of Q, especially in the lower ranges of Q (0 < Q < l). 
This sensitivity is illustrated in figure 3 in which A cr is plotted as a func- 
tion of Q for various values of R x . The steep slopes of these curves for 
0 < Q < 1 indicate the large influence that the core stiffness has on the flut- 
ter boundary in this range. The range of validity of the results in figures 2 
and 3 is limited by the buckling characteristics of the panel. For example, a 
square panel with no stress in the y-direction ^Ry = 0^ and Q = 0.1 buckles 

if R x = 3.3 (for A = 0) . In this case the curve for Q = 0.1 in figure 2 is 
only valid for R x up to approximately 3*3 (see ref. 3)* Other values of a/b 
and Ry yield other ranges of validity. 

It should be noted that all parameters considered in figure 2 (A cr , Q, 
and R x ) are functions of the panel flexural stiffness D s . As a result, a 



Figure 2.- Effect of transverse shear 
stiffness on flutter boundary for 
flat sandwich panel. C =0; 
j = v = 1. 



Q 


Figure 3»- Variation of flutter param- 
eter with transfer shear stiffness 
for several values of streamwise 
in-plane load parameter, v =1; 

Ry — 0; C = 0. 
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sandwich panel having the same flutter speed as a homogeneous panel provides a 
reduction in weight that is not obvious from figure 2. For example, a typical 
flat honeycomb sandwich panel with cover plates and core of the same material. 


a core density of 3 percent of the material density, — = 10, a length-width 

t s 

ratio a/b of 1, length a of 20, core stiffness parameter Q of 0 . 0521 , 
and cover plate thickness t s of 0.044l inch flutters at the same dynamic 


pressure q as a homogeneous panel of the same material and a thickness of 
0.3^2 inch. The weight of the sandwich panel, however, is only 29 percent of 
the weight of the homogeneous panel. 


* 


Curved Sandwich Panel 


Flutter boundaries for a curved sandwich panel are obtained from equa- 
tion (6) . The effect of curvature on A cr is given by the expression 


AA _ (2J + ±) 2 x k n [2j 2 (j + 1 ) 2 v 2 + (2j 2 + 2j + i )v g 

— O ~ 

8j(j + 1) K in 2 


which is the same correction factor derived in reference 2. The curvature 
effects are entirely contained in this expression and are completely separated 
fiom the effects of the core stiffness parameter Q. For this reason, the 
effects of curvature can be thought of in the form of a correction term which 
can be added to the value of A cr for a flat sandwich panel, as is done in ref- 
erence 2 for homogeneous panels. 


One important characteristic of curved panels is that, whereas a flat panel 
has a minimum flutter boundary for j = n = 1 for the range of considered 

herein, a curved panel can have a minimum boundary with other combinations of 
j and n. As a result, various combinations of j and n must be examined 
by trial and error in order to determine the combination that results in the 
minimum flutter boundary. A more detailed discussion on the flutter of curved 
panels is given in reference 2. Although the discussion in reference 2 is con- 
cerned specifically with homogeneous panels, the essential features also apply 
to sandwich panels. 

Figure ^ is a plot of the flutter boundaries for curved sandwich panels 
for various values of the core stiffness parameter Q and for specific values 

- = l\ . The bound- 

* ) 

aries have been minimized (by trial and error) with respect to j and n and 
are, therefore, the lowest boundaries for the range of R x considered. The 

flutter boundary is greatly affected by the core stiffness in much the same 
manner as it was for the flat panel. Flutter boundaries for curved sandwich 

panels = 0.1 and ^ = lj are shown in figure 5 for various values of C. 


of the curvature parameter and length-width ratio (C = 100; 
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The boundaries have again been minimized with respect to j and n. It can be 
seen from this figure that curvature tends to stiffen the panel and raise the 
flutter boundary. 



R y 

Figure 4.- Effect of transverse shear 
stiffness on flutter boundary for 
curved sandwich panel. C = 100; 

* =1. 

b 


t 



Figure 5.- Effect of curvature on 
flutter boundary for sandwich 

panel. — = 1; Q = 0.1. 
b 


CONCLUDING REMARKS 


A small-deflection theory is used in conjunction with a two-mode Galerkin 
solution to analyze approximately the flutter of sandwich panels, both flat and 
slightly curved. Closed-form expressions are obtained from which flutter bound- 
aries can be determined. For flat sandwich panels the flutter boundaries are 
determined in a straightforward manner, but for curved sandwich panels the 
determination of boundaries requires a trial and error minimization process. 

For a given application and flutter speed a sandwich panel can provide a 
substantial reduction in weight as compared with a homogeneous panel. In gen- ^ 

eral, the effect of core transverse shear stiffness on the flutter character- 
istics of sandwich panels is found to be significant. Curvature effects are 
found to be completely separated from core stiffness effects and within the f 

limitations of the present analysis can be treated in the same way as curvature 
of homogeneous panels. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., November 22 , 1963. 
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